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We summarize our recent results for the (T, /x)-phase diagram of Nj = 3 QCD. 
Oh, Finding a strong variation of the critical endpoint fj, c with the quark mass m, we 

point out that an endpoint within reach of experiment requires fine-tuned quark 
' masses. We further discuss the strategy and first data towards extending our 

results to the physical Nf = 2 + 1 theory. 

X 

1. Introduction 



Currently, several approaches are used to side-step the sign problem of 
lattice QCD in attempts to compute the low density QCD phase diagram. 
These are based on multi-parameter reweighting 1 , Taylor expansions 2 and 
analytic continuation from simulations at imaginary fj, 3 , for which there is 
no sign problem 4 . While reweighting is exponentially difficult on large vol- 
umes, calculation of Taylor coefficients becomes increasingly complicated 
with growing order. By contrast, a simulation at imaginary /i = is 
no more difficult than one at fj, = 0. Moreover, it is the only approach 
for which two simulation parameters are varied, resulting in uncorrelated 
statistical errors for different values of /x/. Monte Carlo results for imag- 
inary \i thus involve no approximations. Within the circle of analyticity 
around /i = 0, they can be fitted by truncated Taylor series in (/ii/T) 2 and 
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analytically continued to real fi. While truncating the series introduces an 
approximation, it is well controlled because it can be compared to the full 
result at imaginary /x. 

The task for lattice simulations is to first map out the phase diagram 
in the (m Ul d, fn s ) quark mass plane at [i = 0, i.e. to find the critical line 
rn c s {m u>l i) of second order transitions, separating the first order region from 
the crossover region. In a second step, we wish to know the surface em- 
anating from this line in the /i / direction 6 . So far, even at fx = 
we know only one point of the critical line, corresponding to degenerate 
quark masses, m c (/x = 0),Nf = 3, and only for one coarse lattice spacing, 
N t = 4. In a recent paper we have computed the /i ^ phase diagram 
for the Nf = 3 theory and determined the quark mass dependence of the 
critical endpoint 5 . Here we summarize these results and extend them to 
the physical case of Nf = 2 + 1. The general properties of QCD with com- 
plex chemical potential as well as analyticity properties of observables are 
discussed in the literature 3,5 . 



2. The pseudo-critical line for Nf — 3 

For Nf — 3, we computed the pseudo-critical gauge couplings at the phase 
boundary, and their dependence on imaginary chemical potential as well as 
the quark mass. The data are shown in Figure 1 for different lattice sizes. 
They are well fitted by a polynomial of degree one in the (am) and two in 
(afj,i) 2 . We have explicitly checked that additional polynomial terms are 
statistically insignificant. The figure shows that, to good approximation, 
the leading /ij-tcrm suffices over the interval considered, but our accuracy 
is good enough to also determine the next-to- leading coefficient. Analytic 
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Figure 2. Left: Phase diagram for different quark masses, dotted lines denote crossover. 
The bold curve T*(p?) represents the line of endpoints. Right: Lines of critical quark 
mass separating first order and crossover regions. 



continuation and conversion to physical units yields 

T c ( M = 0,m c (0)) - 1 + 1 - y4 ^ 
+0.602(9) 



/i 



7rT c (0, m< 

2 

+ 0.23(9) 



(1) 



,7rT c (0,m c ); ■ \7iT c (0,m c ), 

Note that the natural expansion parameters are m/ (7rT), ^/(ttT) 2 , leading 
to Taylor coefficients of order one. This is understood in terms of the finite 
temperature Matsubara modes ~ 7rT, which set the energy scale. Since 
the strange quark mass and its difference with m c (0) are small compared 
to 7rT c , the critical temperature is fairly insensitive to quark masses in the 
light quark regime. It is then sensible to compare our Nf = 3 results with 
the Nf = 2 + 1 curve of Fodor and Katz 1 (Figure 1 right), and excellent 
agreement is found. Note that T c (p) decreases only at the percent level 
over the whole range of baryon chemical potentials considered here, which 
thus has a very weak effect. This is a natural consequence of the smallness 
of the effective expansion parameter with Taylor coefficients of order one. 



3. The critical endpoint and its quark mass dependence 

In the three flavor theory, we have a 3d parameter space, (T,fj,,m). The 
pseudo-critical temperature, T c (/x, m), represents a surface in this space. 
On this surface there is a line of critical points, T*(/i 2 ) = T c (/i 2 , m c (fi)), 
separating first order transitions from crossover. Projections of this situ- 
ation to the (T, /j, 2 ), (to, /i 2 ) planes are shown in Figure 2. As the quark 
mass is increased from m c (0), the critical endpoint moves to a finite fi c , 
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while a decrease moves it to imaginary /x c . The change of m c with /x, or 
vice versa, is best seen in the (to, /i 2 )-plane, where it forms a critical line 
separating regions of first order transitions and crossover. We determine 
this line around fi = by measuring the Binder cumulant of ipip. In the 
infinite volume limit it assumes a universal value at a 3d Ising critical point, 

x {(Hib) A ) 

5 4 (to c , Mc )= ,7 2 -> 1-604, K^co, (2) 
(({#)) 

while it is smaller or larger in first order and crossover regions, respectively. 
Data for various quark masses, as well as its universal finite size scaling, 
are shown in Figure 3. Fitting the data to a Taylor expansion about the 
critical point, B±(am, afx) = 1.604 + B (am — am c (0) + A(a/j,) 2 ) + . . ., the 
parameter A represents the first coefficient of the series 



TO c (^0)^ 1 +0 - 84(36) (^) 2 + -- ^ 

Like for T c (/i), the leading Taylor coefficient is of order one. This im- 
plies that the change in to c (/^) is very weak up to chemical potentials of 
fiB = 3/x ~ 500 MeV. Conversely, /i c (m) is extremely sensitive to the pre- 
cise value of the quark mass. We believe that this statement holds for all 
discretizations, because the /z-independent multiplicative mass renormal- 
ization drops out of the ratio in (3), which should have only additive 0(a 2 ) 
corrections. This means that for growing quark masses to > to c (0), the 
critical point ^( c (to) will rapidly shoot to very large values. Hence a small 
fj, c requires a fine-tuned quark mass. 
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Figure 4. Measured critical line separating l.O. and crossover for Nf = 2 + 1 at fj, = 
(diamonds), and finite ^t-endpoints of Fodor and Katz 1 (triangles). 



4. Towards the Nf = 2 + 1 phase diagram 

In the case of physical interest, we have a 4d parameter space 
(T, fi, m u ^d, Tn s ). The critical surface can be determined in a two-step pro- 
cedure. First, we need to complete the (m Ul d, m s ) phase diagram at /i = 0, 
i.e. find the line mg(m tt) d,/i = 0). The results of this step are shown in 
Figure 4. As the light quark mass is decreased, strongly non-linear behav- 
ior is encountered, and linear extrapolations from the Nf = 3 result to the 
physical point are clearly ruled out. In a second step in progress, we com- 
pute the change of this line with finite fii, leading to the critical surface 
iTi'siiTi'u.d, m)- K is intriguing to consider the simulation point of Fodor and 
Katz 1 for physical quark mass values (FK2) in this diagram. This point 
is only slightly displaced from our critical line m c s (m u ^d, [i — 0), but corre- 
sponds to n c ~ 120 MeV. This is consistent with our observation of a rapid 
change of /i c for small variations of the quark masses. Since quark masses 
are very cut-off sensitive, one thus has to expect large shifts of the critical 
point in physical units as the continuum limit is approached. 
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